3D Modulated Spin Liquid model applied to URu 2 Si 2 
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We have developed a 3D version for the Modulated Spin Liquid in a body-centered tetragonal lattice structure 
to describe the hidden order observed in URu2Si2 at To ~ 17.5 K. This second order transition is well described 
by our model confirming our earlier hypothesis. The symmetry of the modulation is minimized for Q = 
(1, 1, 1). We assume a linear variation of the interaction parameters with the lattice spacing and our results 
show good agreement with uniaxial and pressure experiments. 



The fascinating hidden order phase observed in the heavy 
fermion URti2Si2 below To = 17.5 K 0~1|2) is in the center of 
great discussion concerning the origin of the mechanism for 
this second order transition. From thermodynamics proper- 
ties IS a huge entropy quench is observed and this cannot be 
explained by a conventional antiferromagnetic transition due 
to the small value of the magnetic moment in this phase J4| . 

The phase diagram obtained for the URu2Si2 is very inter- 
esting [5-8 1 . In ambient pressure, the system undergoes a sec- 
ond order transition at To = 17.5 K to a phase known as hid- 
den order (HO), with a small magnetic moment /i s» 0.02 hb- 
At very small temperature, a superconducting phase is found 
below T w 1.5 K. Applying hydrostatic pressure or uniaxial 
stress, the system turns into an antiferromagnetic (AF) phase 
with a magnetic moment [i w 0.4 fiB for P x « 0.5 GPa or 
er° « 0.33 GPa, respectively. Bakker et al., [9] showed that 
To increases linearly with the uniaxial stress applied along the 
a axis, and it decreases when the uniaxial stress is applied 
along the c axis. Elastic neutron scattering measurements 
with uniaxial stress applied along [1,0,0], [1,1,0] and [0,0,1] 
directions show a very sensitive variation of the ordered mo- 
ment ifTOl . In-plane stress increases /i, while a perpendicular 
stress does not change the ordered moment at all. The tuning 
parameter to the HO-AF transition seems to be the in-plane 
lattice constant a as shown by Bourdarot et al. JSj . Inelastic 
neutron scattering (INS) measurements show the formation 
of two gaps structure in URu2Si2: one at the incommensu- 
rate wave-vector (1.4, 0, 0) and another at the commensurate 
(1, 0, 0). The second gap is a good candidate for a signature 
of the HO, since the commensurate peak disappears in the AF 
phase HE EHm. 

The various theories that have been presented along so far 
to explain the HO paradigm can be separated into sets of itin- 
erant and localized models. Among the localized models we 
can cite the multipolar models |fT3lFT71 . and among the itin- 
erant models there are different conjecture, indicating that the 
HO-AF transition is either driven by spin density wave lfT8T - 
[20) , by hybridization ET]|22) or by orbital AF fl23]. The great 
advantage of our model is the ability to integrate in a natural 
way the HO and AF already in a realistic localized treatment, 
and to leave open the possibility to include charge fluctuation 



effects later on. 

The Modulated Spin Liquid (MSL) in a 2D square lattice 
was developed in our previous work E4ll in order to explain 
the hidden order phase in URu 2 Si2- The 2D version of MSL 
model provides a simple scenario for the URu 2 Si2, with the 
HO phase resulting from a quantum phase transition where 
the local magnetic moments of the AF phase melt, restoring 
the time reversal symmetry, and preserving the lattice break- 
ing symmetry. One general concept introduced by the MSL 
model is that the modulation of the SL reflects the structure 
of the magnetically ordered phase as shown in FIG.[TJc). The 
transport measurements indicate a continuous transition from 
HO to AF phase. The MSL is in this way the Resonant Va- 
lence Bond relative of the AF phase. 

In the present work, we develop an extension for the MSL 
in a realistic 3D body-centered-tetragonal lattice (BCT-lattice) 
model. Our spin liquid (SL) framework in this 3D system is 
shown to be both realistic and experimentally motivated to ex- 
plain the onset of the HO phase in URu 2 Si2. The BCT-lattice 
structure of URu2Si2 is depicted in FIG.[T] For convenience, 
we use here the tetragonal basis (a, b, c) which contains two 
U-atoms per unit cell of the tetragonal-lattice (T-lattice). Ex- 
perimentally, the lattice parameters are |a| = |b| = 4.124 
A and |c| = 9.5817 A in the HO phase, with less than 1% 
of variation when the system is warmed up to room tempera- 
ture m. 

We start with a Heisenberg model Hamiltonian using the 
standard fermionic representation of quantum spins 1/2, 

H = ^2 Jrr'Xro-XRo-'Xr'ct'Xrv , (1) 

<R,R'),o-(t' 

where the fermion annihilation (creation) operators Xna sat " 
isfy the local constraints Y^a=t 4. Xr ct XRo- = 1- For simplic- 
ity, we consider only the magnetic interaction between near- 
est neighbor sites R and R', as depicted in FIG. [TJb). We 
assume this nearest neighbor interaction to be ferromagnetic 
(^rr' = ^iiitra < ) inside each layer and antiferromag- 
netic ( Jrr' = Jinter > 0) in between adjacent layers. This is 
the simplest and the most natural interaction which can repro- 
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FIG. 1. (a) The crystal structure of URu2Si2. The magnetism 
emerges from the 5/ electrons of U atoms. (b) The BCT-lattice Bra- 
vais structure of the uranium atoms. The intra-layer magnetic cou- 
pling is ferromagnetic, Jj ntra < 0. The inter-layer coupling is antifer- 
romagnetic, Jmter > 0. (c) The structure of the magnetic long range 
order characterizing URu2Si2- 



duce the magnetically ordered phase obtained experimentally 
at high pressure Q: an intra-layer ferromagnetism together 
with inter-layer antiferromagnetism (see FIG.[TJc)). 

Generalizing the procedure of Ref. |24l the Heisenberg 
Hamiltonian ([TJ is decoupled for each bond RR' using appro- 
priated Hubbard-Stratonovich transformations. We find the 
following Lagrangian 

A) = ^ X Ro- I + Ar + (T ^ TOR+z j XRct 
Rer \ z / 

~E Ar + E E kR'xLxRV + C.C.] 

R na (Rei n ,R'ei n±1 ) 



2|v?rr'| 2 _ \ ' ™RWr/ 
(R£i„.R'£L„±i) Jmter (R,R)£L„ mWd 



+E E 



(2) 



where, L n denotes the layer n, and the sum over z refers to 
the nearest neighbors within the same layer. In the follow- 
ing, the Hubbard-Stratonovich fields will be replaced by their 
constant, self-consistent, mean-field expressions, i^rr' = 

-Winter E -<Xr ct XRv) and 771r = Jntra Ect^xLtXR^)- 

Note that this magnetic, intra-layer only, decoupling chan- 
nel leads to a degenerate mean-field system, with each layer 
becoming effectively ferromagnetic, but with an easy axis 
completely decoupled from the other layers. This degener- 
acy does not distinguish an artificially fully ferromagnetic or- 
der from the expected AF order depicted by FIG. [T] (c). A 
more general decoupling scheme would consist in splitting 
arbitrarily the inter-layer interaction, Jj nter = Jsl + Jaf, fol- 
lowing closely the procedure used in Ref. [24] the terms with 
Jaf and Jsl are decoupled in the magnetic and SL chan- 
nels respectively. At mean-field level, the degenerescence is 
lifted by the contribution from the inter-layer part of the lo- 
cal Weiss fields, i.e., the contribution originating from Jaf 
terms. Despite an apparent higher complexity, this general- 
ized mean-field problem is formally identical to the one de- 
scribed originally by the Lagrangian Q. Indeed, consider- 



l 



Er 



ing the BCT-lattice coordination numbers, this general decou- 
pling scheme can be derived at mean-field level from the one 
used here by simply mapping J ntra i— > Jj ntra + 2Jaf, and 
Jnta- ^ Jsl = Jnter — Jaf- As we will see next, what 
is remarkable here, with the BCT-lattice, is that the compe- 
tition between the magnetic and the MSL orders is simply 
tunable by changing the ratio (J in tra + 2 J A f)/( Jnter - Jaf), 
which is qualitatively independent from the arbitrary splitting 
if we take J ntel - = Jsl + Jaf- Therefore, in this work, we 
just assume Jaf = and we consider Jntra/ Jnter as the new 
tuning parameter which is phenomenologically associated to 
pressure variations. 

Experimentally, pressure has a direct effect on the ratio be- 
tween the inter-layer magnetic coupling, J nt er, and the intra- 
layer one, Jntra- This is n °t standard for heavy-fermion sys- 
tems, where pressure variations may often change the local 
energy level of the / electrons. This different phenomenolog- 
ical approach is supported here by a strong experimental evi- 
dence: in URu2Si2, pressure favors a magnetic phase. Here, 
of course, the mechanism is not Doniach-like. 

We introduce the Fourier transform of the fields, Xk<x = 

m k = ^ER e ~* kR ™R< = 
-iq ( R+R ' ) 

' V 2 / (^rr/ , where N is 
The phase factor 6* q = q • R 
is introduced in order to fix the origin of the bond lattice at 
real space position Ro = (a + b + c) /4. 

Hereafter we will concentrate our analysis on the following 
mean-field parameters: the uniform SL, $o = P(p,o,a)/VN, 
the modulated SL, <£>q = <Pq/vN, and the Neel staggered 
magnetization AF, Sq^ = tuq^/vN. There is an impor- 
tant difference with the square lattice MSL [(24}. Here, the 
equivalences between wave-vectors Q for the <p fields do not 
refer to the same Brillouin zones as the ones between Qaf 
for the m fields. This is due to the fact that the magnetiza- 
tion fields are defined on the sites, while the SL fields are de- 
fined on the bonds. The symmetry group of the AF phase 
corresponds to a T-lattice, and Q A f is thus defined modulo 
the first Brillouin zone of the T-lattice. We consider that the 
MSL states, similarly to what happens in the AF phase, sat- 
isfy the T-lattice translational symmetries, although different 
Q modulation vectors are defined modulo a larger Brillouin 
zone, characterizing the long range order with T-lattice peri- 
odicity, but with different intra-shell symmetry breaking (see 
FIG. [2). 

For simplicity, we consider at most the case of one kind 
of modulation at a time in the presence of an homogeneous 
solution. If Q is the wave vector associated with the MSL, 
Q AF is the wave vector associated with the AF. We get: 



^ n ^(R6i>,R'€in + l 

the number of lattice sites. 



tor = 5*Q Ajr e^ AF R , tpq = QoVNSq + QciVNS^+Q, where 
5 denotes the Rronecker delta. Note that a completely equiv- 
alent Ansatz can also be made in the direct bond lattice, sim- 
ilarly to what was done earlier in 11241 . namely, (^rr/ = 

<>R,R'+z I <Pq + fQe ^ v " 4 I. 

We introduce the following reduced notation for the SL 
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(a) (1,0,0) 




(0,0,1) 



(b) (0,0,0) 





FIG. 2. Possible SL long range ordered modulations represented 
on the unit cell of the T-lattice. When it is well defined, the as- 
sociated wave-vector is indicated in reduced coordinates, (h,k,l). 
Dotted lines represent (/Jrr/ = ipo + <PQ, and solid lines represent 
Vrr' = <Po — <PQ bond modulation, (a) Here we show the three 
different relevant odd MSL wave-vectors considered, (b) Respec- 
tively the homogeneous SL, an even wave-vector that we exclude 
(see text), and a mixed modulation which may invoke more than one 
wave- vector (not considered in this work). 



modulation wave-vectors: Q = (h,k,l). The parity of the 
modulation can be obtained from the phase factor e l ^( h - k ^ 1 ) = 
±1. The sign + (— ) characterizes wave- vectors with even 
(odd) parity. The only possible MSL characterized by a sin- 
gle modulation wave-vector have odd parity. The definition 
of parity can be extended to the AF wave-vector. We find 
here that the parity of Q AF is odd (see FIG. [TJc)). For sim- 
plification, we consider here, MSL wave-vectors with either 
1 or 3 modulations. Due to the a — > b symmetry, we fi- 
nally need to compare only three types of modulations (see 
FIG. [2]). All these wave-vectors break BCT-lattice symmetry 
and have the periodicity of a T-lattice: Qi = (1,0,0) also 
breaks a rotation symmetry and characterizes an orthorhom- 
bic lattice, Q2 = (0, 0, 1) may break a mirror symmetry, and 
Q3 = (1,1,1) clearly belongs to the T-lattice group. 

The selection between different modulation vectors Q is 
obtained by comparing the corresponding minimized free- 
energies per site, which is given by 



^(Ao,$o,$Q,£q a 
- A + ^— n$n 



k B T 
N 



l + e 



-mi 



J'm 



I* 



k a=± 
2 



Ql 



J i ntra 



S, 



Qa 



(3) 



where the sum over k is taken over the full Brillouin zone and 
the eigenergies are given by 



with 71, k = cos (%0 cos (4r) cos ( 



'^i 2 k + 4|$ Q | 2 7 2 2 k , Q . (4) 
!f) and 7 2 ,k,Q = 



cos Of + cos (!f + cos (tf + SaV For a given 
modulating wave-vector Q, the staggered magnetization, 
Sq^, and the homogeneous and modulated SL parameters, 
$0 an d ®q are obtained directly from the minimization of the 
free-energy function. The free energy is calculated minimiz- 
ing eq. (0), using Powell's method 11231 , with the auxiliary 
equation to fix the number of n f , 1 = ^ £k E Q=± 1+e ^ng • 
FIG. [3] depicts the behavior of the free energy for the three 
different wave vectors Qi, Q2 and Q3. We find, by varying 
the different physical parameters of the model, that the mini- 
mum is always obtained for the wave vector Q3 = (1,1,1). 
These states, which correspond to the space group No. 134 
P^/nnm also appear to be compatible with the crystallo- 
graphic analysis of URu 2 Si2 made by Harima et al. |[T5l . 
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FIG. 3. Free energy of three modulating vectors Qi = (1,0,0), 
Q 2 = (0, 0, 1) and Q 3 = (1, 1, 1) as a function of: (a) T for t" = 
K and a = ao; (b) hopping for T = 2 K and a = ao; and (c) 
1 — a/ do for T = 2 K and t" = K. The parameters used were 
Jintra(ao) = —6.5 K, Jintei (ao) = 37 K and B — 800 K, as described 
in text. 

A n.n.n. hopping t" may be included in order to phe- 
nomenologically take into account some frustration and intra- 
plane spin-liquid contribution. In this case, the intra-layer 
term is the same for momenta k and k + Q, and the extension 
of our model is directly obtained if we perform the change 
fjj h-> + <"$ cos (k a ) cos (fc fc ). 

Motivated by the pressure experiments which are specially 
dedicated to the anisotropy and the uniaxial effects l8l 4T0l . we 
relate our microscopic interaction with pressure. Bourdarot 
et. al. [8] showed that the uniaxial stress is the relevant vari- 
ational parameter to change the behavior of the system. They 
considered that the deformation is in the linear elastic regime. 
In this work, we propose that our parameters Ji ntra and J; ntel 
also vary linearly with the lattice parameter a. The variation 
can be simply written as Ji nt ei-(a) = ./inter (ao) +Bi(l — a /ao) 
and Jintra(a) = Jintra(ao) - B 2 (l - a/an), where a is the 
value of lattice parameter in ambient pressure. 

Here, the fitting parameters chosen to produce a phase di- 
agram in qualitatively good agreement with experiment are: 
■/inter (an) = 37 K is chosen to obtain Tq = 17.5 K and 
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Ji ntra (ao) = —6.5 K is chosen in order to obtain the best ap- 
proximated value for the critical stress 1 — a/ao ~ 1.45 x 
10 -3 O. To have a good agreement with experiment, we 
also choose the linear coefficient of Ji nte r(a), B\, to have the 
same slope of To, as observed experimentally, and we define 
£> 2 = B\ = B for simplicity. Both Ji„ter(a) and Ji n t ra (a) in- 
crease their absolute values when 1 — a/ao increases. Our 
choice of the interaction parameters variation is, of course, 
a simplified view of the experiment: if we apply an uniaxial 
stress, the in-plane lattice parameters become different. In our 
case both in-plane parameters decrease in the same way. 

The resulting phase diagram is shown in FIG. |4] The varia- 
tion onl-a/oo shows very good agreement with the exper- 
imental results. We define T$ q , T$ and 7s p as the critical 
temperatures for the parameters $q, $o and Sq af , respec- 
tively. Increasing 1 — a/ao, the MSL critical temperature T$ p , 
increases linearly until it reaches the AF ordering temperature 
Tsq and then it goes to zero showing a re-entrance behavior 
due the presence of the hopping t" . The homogeneous com- 
ponent T$ , shows a similar variation, although with a bigger 
amplitude for t" different from zero. In our model, $o persists 
for big values of T or 1 — a/ao, but with a small intensity. For 
the sake of simplicity, we define T$ when $ = 0.6 K. The 
Ts Qaf also increases with 1 — a/ao, but it shows two differ- 
ent behaviors: at first, a fast increase when it is inside the SL 
phase and, a linear increase outside the SL. Here the effect of 
hopping is visible: without this effect, Ts Q will present just 
linear slopes. A first order transition can also be obtained for 
Tsq inside the MSL phase when t" is present. 




FIG. 4. The phase diagram for the modulation Q3 as function of 
deformation (1 — a/ao) and T. Inside the surfaces the mean-field 
amplitudes are different from zero. There is a region where the three 
order parameters coexist. The circle (red), square (blue) and triangle 
(green) lines represent the critical temperatures for the parameters 
$q, $0 and 5q AF , respectively. The parameters used are: t" = 
-2.5 K, J mtra M = -6.5 K, J mtcI {a ) = 37 K and B = 800 K. 

In conclusion, we have developed a modulated spin liq- 
uid model in the realistic three-dimensional BCT-lattice. This 
provides a simple scenario for URu 2 Si2, where the hidden or- 
der results from a quantum phase transition with a very un- 
usual behavior: the magnetic moments of the AF phase melt 



at low pressure, restoring the time reversal symmetry, but the 
lattice symmetry breaking is still present. We analyzed how 
this SL melting in a BCT-lattice can lead to different mod- 
ulation wave-vectors, among which Q 3 = (1, 1, 1) is found 
to be the most stable energetically. The theoretical phase di- 
agram reproduces qualitatively well what is observed experi- 
mentally for URu2Si2- We identify a second order transition 
at To fa 17.5 K and a first order transition from the MSL phase 
to the AF phase at low temperature. The linear dependence of 
■/intra and Jj ntel with the variation 1 — a/ao is a key point of our 
study, which is confirmed by experimental results |8|. Our 
results clearly show that the choice of an appropriate mod- 
ulation vector is crucial for the stability of the MSL phase. 
This could be directly checked experimentally by INS mea- 
surements. By comparing all crystallographic directions one 
could find clear evidence for what this preferable modulation 
might be. Raman scattering experiments could also provide 
another independent check of our results since the orientation 
dependence of Raman spectrum could establish if the modu- 
lation is indeed characterized or not by our Q 3 vector. We 
believe that our study is a very good test for a MSL paradigm. 
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